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SUMMARY 


An investigation of the supersonic area rule has heen conducted hy 
rocket-propelled model tests of zero-lift models of a straight wing- 
hody configuration through a range of Mach number from 0.8 to 1.5 and 
Reynolds number from 5 X lo 6 to 12 X lo 6 . The body of the basic con- 
figuration was modified with axial symmetrical indentations according 
to area-rule considerations for low wave drag at design Mach numbers of 
1.10 and 1.4l. The indentations were designed using a short method 
(given in appendix) that greatly simplified the work. Each configura- 
tion had at its respective design Mach number the same distribution of 
cross-sectional area as the basic body alone. The basic body was 
parabolic having a frontal area equal to O .0606 that of the total wing 
plan-form area. The wing had an aspect ratio of taper ratio of 

0 . 59 ^^ 0° of sweep along the 75-P^rcent chord line, and an NACA 65A004.5 
airfoil section in the free-stream direction. 

Both the Mach number 1.10 and l.kl indentations reduced the wave 
drag of the basic configuration at transonic speeds with the greater 
reduction being obtained from the Mach number 1.10 indentation. The 
beneficial effects from the indentations decreased with increasing 
Mach niimber until Mach number 1.5> above which no benefits were obtained 
from either indentation. Both indented configurations had the same drag- 
rise above Mach number I.I 5 . Tlie theoretical wave drags of the configu- 
rations were useful in showing the relative merits of the different 
indentations . 


INTRODUCTION 


The transonic area rule of reference 1 provides a simple and effec- 
tive means for designing high-speed aircraft for low wave drag near the 
speed of sound. Investigations of the transonic area rule (refs. 1 to ll) 
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for several ving-TDody configurations have shown that a substantial 
reduction in the drag rise near Nhich number 1.0 may be obtained by 
modifying or designing an airplane to have a cross-sectional area dis- 
tribution that is conducive to low drag. A convenient type of modifi- 
cation is a fuselage indentation for Mach number 1.0, which effectively 
cancels the exposed wing cross-sectional areas normal to the axis of 
symmetry. However, recent investigations (refs. 2 and 11, for example) 
have shown that the beneficial effects obtained from such indentations 
decrease with increasing Mach number and eventually produce an unfavor- 
able effect on the drag. This effect is particularly acute on straight 
wing-body combinations where the indentations result in consiaeiauxe 
necking down of the fuselage causing high body slopes at the wing-body 
juncture. 

The concept of the transonic area rule has been extended to super- 
sonic speeds (refs. 12 to 15 ) in an attenpt to provide area distribu- 
tions for low wave drag at supersonic speeds as well as transonic speeds. 
This paper presents an investigation of the supersonic area rule for a 
straight wing-body configuration which was optimized with symmetrical 
body indentations (according to ref. 15) for design Mach numbers of 1.10 
and 1.41. The wing had an aspect ratio of 5-04, taper ratio of 0.594, 

0° of sweep along the 75 -percent chord line, and an MCA 65A004.5 airfoil 
section in the free-stream direction. The unmodified fuselage was para- 
bolic having a frontal area equal to O.0606 that of the total wing plan- 
form area. 

The models were flight tested at the langley Pilotless Aircraft 
Research Station at Wallops Island, Va. The tests of the indented con- 
figurations covered a continuous range of Mach n\mber from 0.8 to I .5 
with corresponding Reynolds n\miber from 5 k 10° to 12 X 10°, based on 
wing mean aerodynamic chord. The experimental results are compared with 
the theoretical wave drags of the configurations tested throu^ most of 
the Nfech nxmiber range. 


SYMBOIS 


A cross-sectional area 

a tangential acceleration 

total drag coefficient, based on S^ 
c mean aerodynamic chord of wing 

g acceleration due to gravity 
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L length of body 

M free-stresLin Mach number 

q free-stream dynamic pressure 

R Reynolds number, based on c 

total wing plan-form area 
W weight of model during deceleration 

X station measured from body nose 

7 angle between flight path and horizontal 

<|> angle of roll of configuration with respect to Mach planes. 

At 0^ of roll the Mach planes are perpendicular to the 
wing plane. 


MODEIS 


Details and dimensions of the models tested are given in figure 1 
and tables I to IV. Photographs of the models are presented as 
figure 2. 

The basic configuration, model A, was originally tested for the 
investigation of reference l6 and consisted of a straight wing on a 
parabolic body with two vertical stabilizing fins. The parabolic 
body profile was formed by two parabolas joined at the maximum diam- 
eter (40-percent body station) and had an overall fineness ratio of 
10.0. The straight wing had a total aspect ratio of 5-0^^ taper ratio 
of 0.59^^ 0^ of sweep at the 75 -percent chord line, and an NACA 65AOO4.5 
airfoil in the free-stream direction. The ratio of body frontal area 
to total wing plan-form area was O.O606. 

Configurations B and C were obtained by indenting the parabolic 
body of the basic configuration for design Mach numbers of 1.10 and 
1.4l, respectively. As is stated in reference I5 for radially symmet- 
rical modifications, the area used for the optimum indentation is 
obtained by averaging the frontal projection of wing areas cut by 
Mach planes at all angles of roll (<5) of the Mach planes with respect 
to the configuration. For symmetrical models, only the average areas 
between 0^ and 90^ have to be considered. Since such computations 
require considerable time, a short method has been devised from which 
the average projected wing areas used for indenting the body can be 
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obtained directly. This method is outlined in the appendix and shows 
how average supersonic area distributions could be obtained by substi- 
tuting a series of hoops for the wing. The resiilting area distributions 
for models B and C, at their respective design I^ch numbers, .were the 
same as the normal cross-sectional area distribution of the original 
parabolic body. It was assumed that the cross-sectional area of the 
original body did not change throxi^ the Mach mmiber range considered. 
The two stabilizing fins used on models A, B, and C were neglected in 
determining the area distributions of the configurations because of 
their thin sections and high sweepback angle. Both models B and C were 
one-half the scale used for model A. The identation losed removed about 
14 percent of the volume of the original body shape. 


TESTS AM) MEASUREMENTS 


All the models were tested at the Langley Pilotless Aircraft 
Research Station at Wallops Island, Va. Models A, B, and C were zero- 
lift rocket-propelled models that were accelerated from zero-length 
launchers to supersonic speeds by single-stage 6-inch ABL Deacon rocket 
motors. Model A, the largest of the group, was propelled by a rocket 
motor installed in its body; whereas, models B and C were propelled by 
boosters, consisting of the rocket motors with stabilizing fins (fig. 5)^ 
that separated from the models after burnout. Velocity and trajectory 
data were obtained from the CW Doppler velocimeter and the NACA modified 
SCR 5^4 tracking radar unit, respectively. A sxirvey of atmospheric 
conditions including winds aloft was made by radiosonde measiarements 
from an ascending balloon that was released at the time of each launching. 

The flight tests covered continuous ranges of Mach nimiber varying 
between Mach numbers 0.8 and 1.5. The corresponding Reynolds nimibers, 
based on wing mean aerodynamic chord, are shown in figure 4 to vary from 
9 X 106 to 25 X 106 for model A, and 5 X I06 to 12 X I06 for models B 
and C through the Mach number ranges covered. 

The values of total drag coefficient, based on total wing plan-form 
area, for all the models were obtained during decelerating flight with 
the expression 


Cd = --^(a + g sin 7) 
qigSw 

where a was obtained by differentiating the velocity-time curve from 
the CW Doppler velocimeter. A more complete discussion of the method 
for reducing the data is given in reference I7. 
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Th© error in total drag coefficient Cjj was estimated to be of the 
order of + 0.0007 at supersonic speeds and ±0.001 at transonic speeds. 

The Mach numbers were determined within ±0.01 throughout the test range. 

RESULTS AND DISCUSSION 


The variations of total drag coefficient Cp with ^fe.ch number for 
the basic and Indented wing -body configurations eire shown in figure 5* 

The drags presented for model A and the fins were published as part of 
an earlier investigation in reference l6. The curves shown in this 
figure represent the basic drag data from the flight tests and are not 
intended for comparison due to the relatively different surface rough- 
ness of the models (obtained from changing model scale) and the different 
Reynolds number ranges of the tests. 

A comparison of the variation of drag rise (Cj) - of the 

models through the Mach number range is shown in figure 6. At Mach 
number 1.0, the drag rise of the basic configioration was reduced by 
55 percent with the Nfeich number 1.10 indentation and by 20 percent with 
the Mach number 1.4l indentation. These results are in qualitative 
agreement with the concepts of the transonic area rule as may be seen 
by comparing the drag rises with the normal cross-sectional area distri- 
butions of the models in figure . Figure 6 also shows that the 
beneficial effects from both indentations decreased with increasing Mach 
number until Mach number 1.5* At higher speeds, the drag rises of the 
basic configuration and the two indented configurations were about the 
same. For the design conditions of the indented configurations, model B 
with the Mach number 1.10 indentation had 25 percent less drag rise than 
the original configuration at M = 1.10, while model C with the Mach 
nmber 1 . 14-1 indentation had about the same drag rise as the original 
configToratlon at M = 1.14-1. Both models B and C gave the same drag rise 
above M = I.I 5 . These comparisons indicate that the M = 1.14-1 inden- 
tation did not optimize the configuration at its design Mach number . 

The supersonic area-rule concept of references I 8 and 15 makes 
possible the calculation of the wave drag of slender wing-body combina- 
tions at supersonic speeds. While this theory when taken to the limit 
of M = 1.0 shows that the drag is solely dependent on the normal 
cross-sectional area distribution (transonic area rule of ref. l) , the 
magnitude of the theoretical wave drag near M = 1.0 is in error as is 
shown in references 15 to 15 . At supersonic speeds, reference 15 shows 
that the wave drag may be determined from a series of area distributions 
or equivalent bodies of revolution. Each area distribution of the series 
is obtained from the frontal projection of the area intercepted by 
parallel Mach planes at a given angle of roll of the configuration with 


6 


KACA RM L55C09 


respect to the Mach planes. The drag coefficient of the configuration 
is then determined from the average of the drag coefficients for all the 
equivalent bodies of revolution. It is important to use the series of 
area distributions for determining the wave drags and not the average 
area distributions that were employed in designing the indentations. 

Use of the average areas for the computations would give values that 
greatly underestimate the drag of the configurations. A detailed de- 
scription of this method showing the evaluation of the wave drag by 
Fourier sine series and how to determine the wave drags over a range of 
Mach nijmber is given in references lU and 15 . 

Since the present models were symmetrical, only the areas between 
roll angles of 0° and 90° "to be considered. Some of the area dis- 
tributions for the models tested are shown in figures 7(b) to 7(g) . 
According to the convention used, the Ifech planes were perpendicular to 
the wing plane at 0° roll. For the roll angles at which the Mach line 
slices were alined with the blunt leading edge and trailing edge of the 
wing, discontinuities appear on the slopes of the area distributions, 
limiting the application of the linearized theory to about Nfeich number 1.1 
for the present case. Hl^ values of slope were assumed rather than 
these discontinuities in order to extend the Nfeich number range and 
determine if reasonable agreement could be obtained with the experimental 
drag rises at higher ffech numbers. The Fovirier sine series used for 
calculating the drag of each equivalent body of revolution was evaluated 
to 2k terms and plots of these series indicated that the series were 
convergent. It should be noted, however, that while the series converge 
for the first 2k terms they may diverge when more terms are evaluated. 

In this regard, reference 15 suggests that 48 terms be used before 
establishing the convergence of the series. 

The theoretical wave drag coefficients were computed for the models 
thro ugh a range of Mach number from 1.0 to 1.4 and are presented in 
figure 8. The computed wave drags in figure 8 are in general agreement 
with the observations made from the flight test data and lead to the same 
conclusions, namely that the Mach number 1.10 indentation is better than 
the Mach number 1.4l indentation at transonic and low supersonic speeds 
and that the Mach number 1.4l Indentation offers little or no savings in 
wave drag at the design Mach number. Near Mach number 1.0, the theoret- 
ical variations are not realistic due to the limitations of the theory 
(ref. 15 ) used. 

Comparisons of the theoretical wave drag coefficients and the 
measured drag-rise coefficients for each configuration tested are given 
in figure 9 . The best agreement is shown for model B in both magnitude 
and variation. Although the agreement obtained for models A and B is 
poor, the deviation between the theory and test values is not greater 
than 20 percent in the Mach number range from 1.1 to 1.4. The agreement 
obtained in this ^fech number range is better than what was expected in 
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view of the fact that the theory does not apply above M = 1. 1 because 
the leading edge of the wing is blunt. However, from the standpoint of 
the configuration as a whole, the theory was useful in showing the rela- 
tive merits of the different indentations. 


CONCLUSIONS 


The results of an investigation of the supersonic area rule by 
rocket-propelled model tests of zero- lift models of a straight wing- 
body combination having body Indentations for design Mach numbers of 
1.10 and l.Ll, tested through a range of Mach number from 0.8 to 1.5 
indicate the following conclusions: 

1. Both the Mach number 1.10 and 1.4l indentations reduced the wave 
drag of the basic configuration at transonic speeds with the greater 
reduction being obtained from the Mach number 1.10 indentation. Both 
indented configurations had the same drag rise above Mach n\miber 1.15* 

2. The beneficial effects from the indentations decreased with 
increasing Nkich number until Mach number 1.5^ above which no benefits 
were obtained from either the Mach number 1.10 or 1.4l indentations. 

5 . The theoretical wave drags of the configurations were useful in 
showing the relative merits of the indentations. 


Langley Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
langley Field, Va., Febrmry 25, 1955- 
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APPENDIX 

A SHORT METHOD OF DETERMINING THE AVERAGE AREA DISTRIBUTION 
FOR AIRCRAFT AT SUPERSONIC SPEEDS 
^ Maxime A. Faget 


The equivalent average area "body for supersonic speeds is deter- 
mined from slices made at the Mach angle. The area used is the frontal 
projection of the average area cut hy a Mach plane as the aircraft con- 
figuration is rotated about its axis a full 5^0° relative to the Mhch 
plane. 



The basis of the method to be described is to construct an imag- 
inary model of the configuration that represents all rotational positions. 
Th\is, one plane slanted at the tfech angle will cut an area that is 
already averaged for all positions of relative roll between the actual 
configuration and the cutting plane. This imaginary model consists of 
a body of revolution representing the original fiiselage surrounded by a 
number of hoops representing wings, nacelles, tail fins, and so forth, 
external to the fuselage. The axial distribution of cross-sectional 
area of the body of revolution is eq\ial to that of the original fuselage, 
the hoops represent the axial distribution of cross-sectional area 


Area 
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of the original wings ^ nacelle tail fins^ and so forth. The radius of 
each hoop is equal to the mean distance from the airplane center line 
of the portion of wing^ nacelle, tail fin, and so forth, that the hoop 
represents. Thus, a wing may he divided into 10 spanwise portions. 

Each portion would then be converted into a hoop equal to it in cross- 
sectional area distribution and with a radius equal to the mean distance 
of the panel from the center line. 



The hoops thus obtained, which are called parent hoops, are then 
subdivided into smaller elementary hoops in order that the variation of 
cross-sectional area distribution may be introduced into the computa- 
tional procedure. Each family of hoops thus obtained has the same 
radius as its parent hoop. For convenience, these elementary hoops are 
made to have equal frontal area and rectangular sections, but various 
chord lengths. 
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Thus, these new hoops may be described by leading-edge station, 
trailing-edge station, and hoop radius. It is not necessary to 
designate the hoop thickness as this is fixed by the fact that all hoops 
are of equal frontal area. Thus, the thickness is inversely proportional 
to hoop radius (hoops in the same family have the same thickness) . With 
hoops of this type geometry, a simple relationship exists between the 
projected frontal area of the hoop which is cut by the Mach plane and 
the relative positions of the Mach plane and the hoop, so that rapid 
computation of the areas may be carried out. 

fFhp following notations and symbols are used to describe the ele- 
mentary hoop geometry and the Mach plane cutting the hoop: 

total frontal area of the hoop 

frontal projection of hoop area cut by plane 

total circumference of hoop 

partial circumference of hoop leading edge ahead of Mach 
plane 

partial circumference of hoop trailing edge ahead of Mach 
plane 

radius of the hoop 


A* 

A 

C 

^le 

Cte 

R 


R' 


R 

tan p 


L 

T 

X 


hoop leading-edge station 
hoop trailing-edge station 

distance from the nose to the intersection of the Mach 
plane and the center line 


Mach angle 
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From the geometry in the preceding sketch, the following may he 
deduced: 


A _ 

A* C ' C 


Cle _ 
C 


^TE _ 
C 


cos' 


\( (L-X) tan p' 


R 


cos 


180 

_T ^(T-X)tan 

180 


Note: 

„han ilcaaiLtt > 1 , let cos-1 , q 

„hen K < . 1 , let coe-l( < i rJ H .') = 180 

R \ R j 


( 1 ) 

( 2 ) 

( 5 ) 


Actual construction of the hoops is unnecessary for computation 
of the effective area distribution. The wing is divided into spanwise 
panels and the normal cross-sectional -area variation of these panels is 
then determined and plotted against model station. 
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Each individual area plot is then further broken down into a 
number of rectangles, each rectangle being equal in frontal area. The 
leading and trailing edges of the rectangles are chosen to best repre- 
sent the area progression of the particular area plot. These rectan- 
gles represent the elemental hoops upon which the computations are done 

The computational procedures may be carried out on a simple compu- 
tation form as is shown below. The leading -edge and trailing-edge 
positions and the hoop radii are first reduced to dimensionless param- 
eters. A sample conputation showing this procedure is given for two 
families of hoops for a M = 1.82 (tan p, = 0.66). 



Station, in. 


Panel "G" 
R = l6" 



Station, in. 
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The area cut at any station X is next determined by first 

computing values of ~ ^ and T - X and then determining Ct^/C 

R' R' 

and Cte/C from a plot of these functions (fig. lO) or from equa- 
tions (2) and ( 5 ) . The computation for this is illustrated in the 
following table for station X = 25* 



The average area cut by the Mach plane passing throu^ station 25 
is then determined by summation of all values of Cte/C and Cte/C. 
The difference between these summations multiplied by the frontal area 
of the elemental hoops ( 0.25 sq in. for this case) is the projected 
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frontal area cut ty the Mach plane. This may he expressed as follows. 


A = A* 




Cte 

c 


For the design of indented configurations, according to the area- 
rule concept, the indentations are obtained by subtracting the summed 
projected area of the wing from the body cross-sectional area at 
corresponding stations. 
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TABLE I 

COORDIMTES OF MCA 65A004.5 AIRFOIL 


Station^ 

Ordinate, 

percent chord 

percent chord 

0 

0 

.5 

• 5i^9 

• 75 

.421+ 

1.25 

.540 

2-5 

.758 

5.0 

.986 

7-5 

1.194 

10 

1.568 

15 

1.646 

20 

1.855 

25 

2.014 

50 

2.151 

55 

2.208 

1+0. 

2.246 

1^5 

2.245 

50 

2.196 

55 

2.099 

6o 

1.957 

65 

1.780 

70 

1.572 

75 

1.558 

80 

1.084 

85 

.818 

90 

.549 

95 

.280 

100 

.010 

L.E. radius: O.I5O percent chord 

T.E. radius: O.OII5 percent chord 
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TABLE II 


COORDIMTES OF PARABOLIC BODY^ 


[stations measured from "body nosej 


Station, 

in. 

Ordinate, 

in. 

0 

0 

1 

.245 

2 

.kQl 


.925 

6 

1.527 

10 

2.019 

Ik 

2.558 

18 

2.942 

22 

5.175 

26 

5.250 

50 

5.255 

3 k 

5.181 

58 

5.095 

k 2 

2.975 

k 6 

2.820 

50 

2.651 

3 k 

2.407 

58 

2.149 

62 

1.857 

65 

1.615 


iBased on length of the 65 -inch 

hody . 
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TABLE III 

COORDIMTES OF BODY WITH MACH NUMBER 1.10 INDENTATION 
[^Stations measiared from Body nosej 


Station, 

in. 

Ordinate, 

in. 

(a) 

(a) 

50 

5.255 

52 

5.160 


2.920 

56 

2.650 

58 

2.575 

40 

2.185 

42 

2.095 

44 

2.108 

46 

2.185 

48 

2.272 

50 

2.548 

52 

2.402 

54 

2.575 

56 

2.285 

58 

2.149 

60 

2.007 

62 

1.857 

64 

1.698 

65 

1.615 


^Coordinates Between stations 0 
and 50 are identical to those of 
taBle II. 
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TABLE IV 

COORDINATES OF BODY WITH MACH NUMBER 1 . 4 l INDENTATION 
[^Stations measured from "body nose^ 


Station, 

in. 

Ordinate, 

in. 

(a) 

(a) 

22 

5.175 

2 k 

5.190 

26 

5.175 

28 

5 .i 4 o 

50 

5.080 

52 

2.985 

54 

2.870 

56 

2.742 

58 

2.642 

40 

2.560 

42 

2.482 

44 

2.422 

46 

2.568 

48 

2.525 

50 

2.288 

52 

2.252 

54 

2.220 

56 

2.152 

58 

2.042 

60 

1.952 

62 

1.790 

64 

1.650 

65 

1.610 


^•Coordinates "between stations 0 
and 22 are identical to those of 
table II. 



Model Characteristics ro 

Wing aspect ratio 5.0I4. 

Wing taper ratio 0,59i4. 

Wing mean aerodynamic chord, ft 2.575 

Free stream airfoil, NACA 65AOOI4..5 .. .Table I 

Total wing planform area, sq ft I5.208 

Body fineness ratio 10.0 

Body frontal area, sq ft 0.922 

Body coordinates. Model A Table II 


^5.5H 



A -A 

Typical fin section 



(a) Basic wing-body combination (ref. I 6 ). Model A. 


Figure 1.- 


Details and dimensions of models. 


All dimensions are in inches 


MCA RM L55C09 



Model Characteristics 

Wing aspect ratio 3 •Oh 

Wing taper ratio.'. 0.59^4- 

Wing mean aerodynamic chord, ft I.I87 

Free stream airfoil, NACA 65AOOI4..5 .. .Table I 

Total wing planform area, sq ft 

Body coordinates: Model B Table III 

Model C Table IV 
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A-A 

Typical fin section 



(b) Configurations with body indentation for Mach number 1.10 and l.^t-l. 

Model B and Model C. 


ro 

UJ 


Figure 1.- Concluded 
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(a) Basic configuration (ref. l6). Model A. L- 721 I 4 . 7 .I 



Figure 2.- Photographs of models. 




> 



(c) Conf igioration Indented for M = 1.4l. 


Model C. L-81+i|.53,l 


ro 

VJl 


Figure 2.- Concluded 
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L-8i|685.l 

Figure 3»- Photograph of a model and booster on zero-length launcher. 
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Figure k.- Variation of Reynolds number with Mach n-umber. Reynolds number 
is based on the mean aerodynamic chord of the wing. 



.8 .9 1.0 1.1 1.2 1.5 1.4 1.5 


Figure 5*- Variation of total drag coefficient with Mach number for models 

tested. 
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Figure 6.- Comparison of the variation of drag rise coefficient with 

Mach number for models tested. 
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Mach angle 


Roll angle, ^ — ' 




(a) Cross-sectional area distributions. M = 1.0. Models A, and C. 

Figure 7«- Comparisons of area distributions of models tested at various 

Mach numbers. 


o 
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(b) Basic configuration at M = 1.10. Model A. 



(c) Configuration indented for Mach number 1.10 at M = 1.10. Model B. 



(d) Configuration Indented for Mach number 1.4l at M = 1.10. Model C. 


Figure 7*- Continued. 
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(f) Configuration indented for Mach mmber 1.10 at M = 1.4l. 



Model B. 


(g) Configuration indented for Mach number 1.4l at M = 1.4l. Model C. 

Figure 7*- Concluded. 



Figxir© 8.— Comparison of "the variation of the theoretical wave drag 
coefficients with Mach number for models A, B, and C. 
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(a) Basic conflg\iration. 



(b) Configuration indented for M = 1.10. 



(c) Configuration indented for M = l.Ul. 


Figure 9 .- Comparison of the measured drag rise coefficient and the 
theoretical wave drag coefficient for each model tested. 
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Figure 10.- Variation of area parameters for use with computational 

procediire given in appendix. 
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